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FLOWS AND STOCHASTIC TAYLOR SERIES 
IN ITO CALCULUS 

KURUSCH EBRAHIMI-FARD“, SIMON J.A. MALHAM\ 
FREDERIC PATRAS= AND ANKE WIESE** 


Abstract. For stochastic systems driven by continuous semimartingales an 
explicit formula for the logarithm of the Ito flow map is given. A similar 
formula is also obtained for solutions of linear matrix-valued SDEs driven 
by arbitrary semimartingales. The computation relies on the lift to quasi- 
shuffle algebras of formulas involving products of Ito integrals of semimartin¬ 
gales. Whereas the Chen—Strichartz formula computing the logarithm of the 
Stratonovich flow map is classically expanded as a formal sum indexed by 
permutations, the analogous formula in Ito calculus is naturally indexed by 
surjections. This reflects the change of algebraic background involved in the 
transition between the two integration theories. 


1. Introduction 


The setting of our work is classical stochastic calculus, as exposed, e.g., in the 
classical textbooks [20l[2^. Its aim is to obtain an explicit formula for the logarithm 
of the Ito flow map associated to a stochastic differential system, generalizing the 
classical work of Ben Arous [5] on Stratonovich flows and stochastic Taylor series. 

In this introduction, we briefly state the main result of this work. A description 
of the historical background and details on the underlying definitions and objects 
of study are postponed to the next section. 

Let {X^, X ^,..., X^} be scalar continuous semimartingales. We assume, with¬ 
out loss of generality, that Aq = 0 and that their quadratic covariation, or square 
bracket operation, is such that [X'^,X^] = 0 for i ^ j. We consider the general 
stochastic differential system (in what follows the notation /q • • • dX] refers to Ito 
integrals) 


N 


( 1 ) 


U = Yn 


/ v.iYs)dx: 


where U : t are smooth vector fields. In the following, we will identify V) 

with the partial differential operator Vidy := J2j=i For i = N + 1,..., 2N, 

let A* denote the quadratic variation of X’’~^, that is A* = A*"-^], and 

define the second order differential operator V) by 


( 2 ) 



E 


j, k=l 
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By analogy to the Stratonovich stochastic system driven by Wiener processes in 
[I], we define the stochastic partial differential operator S as 


2N 

i=l 


Since the differential operators Vi are time-homogeneous, we also have dSt := 
For n > 1, let / S'” denote the n-times repeated integral of S and set 
/ SO = Id, so that 


/ 


gn 



2N 

Vjn 

il Jn = l 


dX^ . . .dX/" 


The Ito-Taylor series expansion for the ffowmap ipt '■ Yq ^ Yt corresponding to the 
stochastic system is given by 

<pt = j2 loi. 

n>0'^ 


We set S := Y.n>o / S” = Id -I- / SdS, which describes the action of the ffowmap 
on smooth functions. The central aim of this work is to calculate logS in the Ito 
framework. 

The computation of the logarithm of this action may be considered as the sto¬ 
chastic analog of a well-known problem in the theory of classical differential equa¬ 
tions (motivated, e.g., by numerical considerations, and referred often to as the 
continuous Baker-Campbell-Hausdorff problem). And it turns out that, provided 
one uses Stratonovich integrals, the solution for stochastic differential equations is 
essentially the same as in the classical case. We refer the reader to the next section 
for more details. 

Displaying the first few terms of log S in the Ito framework may give an idea of 
the complexity of its expression. We adopt the notation = [X*,X-^], and for 
a repeated Ito integral /ji,..., := / ••• / dX^ ■ ■ ■ dX^. Then the first three terms 

are as follows 


logS= + 

i i,j 








^kij 


i,j]k 


+ Ijli,k]) 


I 

6 


{likj + Ijki + Ii[j,k]) + -^{lkji + I[j,k]i + ^kli, 



To write the general expression for this expansion, let us introduce some notation. 
For / a surjection from the set [n] n} to the set [fc] (written / G Sj^ f,), 

we set 


d{f) ■■= |{* < n,f{i) > f{i + l)}\. 

The set of surjections / from [n] to [fc] such that Vz < k, |/“^{z}| < 2 is written 
fc For a sequence J = (ji,... ,jn) of elements of [2X] and A = Ai\[ - ■ ■\[Ak = 
[n] an ordered partition of the set [n] into disjoint subsets, we write for the 
iterated integral / ■ ■ ■ / ■ • • dX^*‘ , where, for Ai = {m,..., a;}, X/* stands 
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for the iterated quadratic covariation ,..., := , [X^‘^ 2 , X^'^i\[. 

For / as above, we set A{f) := U • • • U f~^{k) and 

2N 

sj:= Y. 

ii . ., 2 „ = 1 

Our main result reads 


Theorem 1.1. IFe have: 


iog(s) = E E E 

ra>0 n>k>l 


i-YYl 

n 



Sf- 


This statement follows from Theorem [621 The restriction of the indexing set to 
Sj)^ fc follows from the fact that, for continuous semimartingales, iterated brackets 
such as ^Xi,Xj\,X}^] vanish. When allowing semimartingales with jumps (the 
article will develop the combinatorial theory of iterated integrals of semimartingales 
in this more general setting), this restriction disappears. 

For example, with / the surjection from [3] to [2], defined by 

/(I) =2, /(2) = 1, /(3) = 2, 

we obtain d{f) = 1, and 

(_ 1 ) 4 /) (- 1 ) / 2 \-^ 1 

n \d{f)) 3 \l) 6 

and Sf = ViVjVkIj[ik], as expected from the low order direct computation given 

i,j,k 

previously. 

The rest of the paper develops the formalism necessary to prove the above the¬ 
orem. Several of the tools that enter our approach are of general interest, and 
allow to handle the algebraic structures of iterated ltd integrals. We also show how 
the same ideas can be applied to the study of linear stochastic matrix differential 
equations driven by arbitrary semimartingales fTheorem 17.31) . 

Acknowledgements: The first author is supported by a Ramon y Cajal re¬ 
search grant from the Spanish government. The third author acknowledges support 
from the grant ANR-12-BS01-0017, Combinatoire Algebrique, Resurgence, Monies 
et Applications. 


2. The Strichartz formula 

Let us recall first the historical as well as technical background of Theorem ll.il 
Its knowledge will help to enlighten our forthcoming constructions, and make clear 
to what extend ltd calculus differs from the Riemann or Stratonovich calculus. 

From the seminal 1957 work by K.T. Chen on the algebraic structures underlying 
products of iterated integrals [B] followed the existence of an exponential solution 
of the classical nonautonomous initial value problem 

Y{t)=F{t,Yit)), Y{0)=Yo, 

where F{t) : —>■ is a vector field depending continuously on time. 

Over the decades following Chen’s work, the explicit formula for this exponential 
solution was obtained independently by several authors. Mielnik and Plebanski [53] 
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as well as Strichartz [29] calculated the function il{t) such that Y(t) = exp(f2(t))yb- 
Using the fact that Stratonovich integrals obey Chen’s rules of calculus for iterated 
integrals, Ben Arous showed soon after the work of Strichartz, that the application 
domain of this formula extends to stochastic realm mm- 

The explicit expression of the series 17 (t) is rather intricate. Its most classical 
formulation involves permutations and iterated integrals of iterated Lie brackets 

= E E - 

n>0a&Sn ^ 

where the bracket of vector fields follows from their interpretation as differential 
operators, and the integration domain is the n-dimensional simplex 






t/. ' 

■)/ ^[0,t] 


[A'(Scr(l) ), F{^S(j(^2) )]? * * * ' * ' dSn, 


^fo.d ■ • ■ ’ 0 ^ ■Si < • • • < s™ < t}- 

In the formula for 17(7), S'„ denotes the set of permutations of the set [n] := 
n}. The quantity d{a) (already introduced in the more general case of 
surjections) is called the number of descents of the permutation a € Sn, that is 
the number of positions in the permutation (cr(l),..., cr{n)), where a{i) > a{i + 1), 
for i = 1,... ,n — 1. The formula is known in the literature as Strichartz, Chen- 
Strichartz or continuous Baker-Campbell-Hausdorff formula. 

It can also be stated by replacing iterated Lie brackets by standard operator 
products. The formula is then essentially the same except for the scalar coeffi¬ 
cient, which becomes -the same kind of scalar coefficient as 

those appearing in Theorem ll.il We refer to [29] for details on the analytic back¬ 
ground, and to [minj for the underlying combinatorics of free Lie algebras and 
Lie idempotents. 

The derivation of the formula for 17(7) relies on a precise understanding of the 
calculus of iterated integrals, which in turn is based on integration by parts. For 
two scalar valued indefinite integrals F(7) := /(s)ds and G{t) := g(s)ds, recall 
that ^ 

f(s)G(s)ds+ [ F{s)g{s)ds. 

Jo 

For general iterated integrals 

Fn{t) := [ /l(si)/2(S2)---/n(Sn)dsi---ds„, 

the above product generalizes to Chen’s shuffle product formula 


F(7) . G(7) = 


(3) F„(7)-T^(7) 



/o-i(l)('Sl) ■ ■ ■ fa (n+m) (Sn-I- m )dsi■• ■ ds„+ 

m-i 


where Shji_m is the set of (n, to)- shuffles, i.e., permutations a on the set [n + to], 
such that (t( 1) < • • • < cr(n) and a{n -I- 1) < • • • < a{n + to). 

The above product is conveniently abstracted into an algebraically defined shuffle 
product on words. Let Y := {?/i, j/ 2 ,...} be an alphabet, and Y* the correspond¬ 
ing free monoid of words w = ■ ■ ■ yin- The vector space K{y), which is freely 

generated by T*, is a commutative algebra for the shuffle product: 


( 4 ) 


Vi- ■ -VpiaVp+i- ■-Vp+q-.= 

1 p+g 

o-eShp,g 
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with Vj G Y, j G {1,... + g}. We define the empty word 1 as unit: 1 mu = 

u m 1 = u for V G Y*. This product is homogenous with respect to the length of 
words and can be defined recursively. Indeed, one can show that: 

( 5 ) Vi---VpUl Vp+I ■ ■ ■ Vp+q =Vi (U2 • ■ • Up m Vp+I ■ ■ ■ Vp+q) 

+ Vp+I (ui ■ • • Up m Up+2 • • • Up+,). 

The shuffle product was axiomatized in the 50’s in the seminal works of Eilenberg- 
MacLane and Schiitzenberger lillH], and has proven to be essential in many fields 
of pure and applied mathematics. In [Hi Chen studied fundamental groups and 
loop spaces. In control theory, Chen’s abstract shuffle product plays a central 
role in Fliess’ work [14], which combines iterated integrals and formal power series 
in noncommutative variables into an algebraic approach to nonlinear functional 
expansions. Reutenauer’s monograph on free Lie algebras m imbeds Chen’s work 
into an abstract Hopf algebra theoretic setting. More recently, Chen’s formalism 
came to prominence in Lyons’ seminal theory of rough paths |22j . 

When dealing with iterated ltd integrals of semimartingales, this machinery of 
shuffle products does not apply any more. One has instead to use the quasi-shuffle 
product [5] . Its definition will be recalled further below. We simply mention for the 
time being that, although it was discovered and investigated much more recently, 
it appears to be as important as the shuffle product, both from a theoretical as 
well as applied point of view. Indeed, it encodes the algebraic structure of discrete 
sums, as the shuffle product encodes the one of integration maps, and appears 
in various domains, e.g., multiple zeta values, Rota-Baxter algebras, and Ecalle’s 
mould calculus. The latter is related to the computation of normal forms in the 
theory of dynamical systems. We refer to [III na [mill] for further historical and 
technical details on quasi-shuffle algebras. 

3. Semimartingales 

Recall that a process X is a semimartingale, if X has a decomposition Xt = 
Xq + Mt + At for t > 0, where Mq = aIq = 0, and where M is a local martingale 
and A is an adapted process that is right-continuous with left limits and has finite 
variation on each finite interval [0,t]. 

It is well-known that the space of semimartingales with multiplication forms 
an associative algebra. The quadratic covariation or square bracket process [-Y, T] 
between two semimartingales X and Y, is defined via their product as follows 

(6) X • y = XoYo + J X_dY + jY.dX+ [X, Fj. 

The quadratic covariation of a process X with itself is known as its quadratic 
variation. Let X, Y and Z be semimartingales, the quadratic covariation satisfies: 

(1) [X,0] =0; 

(2) [x,y] = [y,x]; 

(3) [X,[y,z]] = [[y,y],y]. 

Hence, the square bracket process defines a commutative and associative product 
on the space of semimartingales. We refer to the monographs by Protter [26] and 
Jacod & Shiryaev m for details. 
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For notational convenience, we will write from now on iterated integrals of semi¬ 
martingales as follows: 


XV := / dr, and 


Jx,-Xr,:= JiJ Xi-.-X„_i)_dX„. 


Iterated brackets are denoted by: 

(7) X*Y:=[X,Y], and :=[■•■ [Xi,X 2 ], ••• ,X„]. 


From now on we will assume that all processes are normalized so that Xq = 0. 
Terms such as XqYq, can therefore be ignored in products of stochastic integrals, 
e.g., as in equation 

Let us briefly illustrate the combinatorial nature of iterated Ito integrals, before 
turning to the general, and more abstract picture. Recall first that, for arbitrary 
semimartingales A,B,C,D and X := f A- dB, Y := f C_ dD, we have 


(8) X-Y = J {XC)_ dD -b y {AY)_ + J [AC)_ d[B, D], 

so that, for example, for B := J dB, Y := f C- dD 


B-Y = J{CB)_ dD +J Y_ dB +J C- d[B,D] 

= J CDB + J{J C- dB + J B_ dC+ [B,C])- dD + j C- d[B,D] 
= J BCD + J CBD + J CDB + J[B,C]D + J C[B,D]. 


Whereas the first three terms of the expansion, i.e., / BCD + J CBD + f CDB, 
are those that would also appear in the shuffle product expansion, the last two 
terms arise from the bracket operation. These extra terms are typical outcomes of 
what distinguishes quasi-shufhe and shuffle computations. The algebra underlying 
the quasi-shuffle calculus has been explored in MM- In a nutshell, the bracket 
terms that appear in the above product require the replacement of permutations in 
the calculation of the Chen-Strichartz formula by the larger class of surjections. 


4. Quasi-shuffle algebra 

We recall now the definition and basic properties of quasi-shuffle algebras. In 
spite of the fact that such algebras encode naturally Ito integral calculus, they 
appeared only in a few papers in the context of stochastic integration, see e.g. the 
works by Gaines, and Liu and Li [nisi] and [5]. As the name indicates, quasi¬ 
shuffle algebras are obtained as a deformation of classical shuffle algebras. It is 
generally agreed that the idea can be traced back to the work of P. Cartier on free 
commutative Rota-Baxter algebras [3]. However, it was formalized only recently by 
M. Hoffman [TH]. The link between Hoffman’s ideas and commutative Rota-Baxter 
algebras was explored in m- 

Abstractly, a quasi-shuffle algebra is defined as a commutative algebra B with 
product •, which is equipped with two extra bilinear products denoted f and j, 
(known as “half-shufhes”), such that for x,y,z G B one has 

(9) xty = yix, 

(10) {x»y)\z = x»{ytz) 


EXPONENTIAL LIE SERIES 


7 


( 11 ) {x^y)^z = x^{y^ z + yiz + yz). 


One writes usually 

x\^y := X "[ y + X i y + X • y, 

and calls UJ the quasi-shuffle product. In particular the last axiom then reads 
{x y) z = X {y\^z). Note that when the product • on B is the null product, 
then one recovers the usual axioms for shuffle algebras [51HH] • 

Let us mention that the “deformation” induced by the • product can be under¬ 
stood in terms of - weight-one - commutative Rota-Baxter algebras. The term 
y^z + y],z + y»z can then be interpreted as the so-called double Rota-Baxter 
product. We refrain however from developing these ideas here, since they are only 
indirectly relevant to our present purposes. The interested reader is referred to the 
survey paper m for an overview of the links between the theories of Rota-Baxter 
algebras, integral calculus and (quasi-)shuffle algebras. 

The standard example of a quasi-shuffle algebra, studied in detail in US], is 
provided by the tensor algebra T{A) := over a commutative algebra 

(A,*). The three products t) • are defined inductively by: a •b := a * b, 
a 'I b := ba, a b := ab, and 


«! ■ 

■ T ' 

' bra 

:= (ai • 

• a„_iLU6i 

■' ’ bm^^n j 

ai ■ 

• dn bi * 

' bm 

:= (ai • 

• a„LU^i • • 

bm — l)bm-, 

Ol • 

••an^bl• 

' bm 

:= (ai • 

• a„_iyJ6i 

■ ' * bm — l){(^ 


where we used the common word notation oi • • ■ a„ for oi ( 8 > • • • ( 8 ) G A®". 
For example, the product of two words of length two gives explicitly 

aia2L!J6i&2 = (ai UJ5162)02 -I- (aia2LU&i)^2 + (aiLU6i)(a2 * 62) 

= 61620102 -I- 016162O2 -I- 61O162O2 -I- 01O26162 -I- 61O1O262 -I- 01610262 
+ 61 (oi * 62)02 + (oi * 61)6202 + Ol(o2 * 61)62 + (oi * 61)0262 
+ 0161(02 * 62) -f 6101(02 * 62) + (oi * 6i)(o2 * 62). 


Recall from Section [5] that the recursive description of shuffle product is comple¬ 
mented by its definition in terms of permutations. Similarly, the above recursive 
definition of the quasi-shuffle product has an explicit presentation in terms of sur¬ 
jections. Concretely, let / be a surjective map from [n] = {!,...,n} to [p\. We 
set: 

/(oi • •-On) := ( Oj)(g)-'-®( 

so that for / from, say, [ 4 ] to [2] given by /(I) = 1, /(2) = 2, /(3) = 1, /( 4 ) = 2, 
we find /(oi • • • 04 ) = (oi * 03 ) ( 8 ) (02 * 04 ). 

Then, we obtain: 

(12) oi • ■ •o„yj6i • •-6^ :=y^/(oi---o„6 i---6 ^), 

/ 

where / runs over all surjections from the set [n-l-m] to the set [fc], for max(n, m) < 
k <m + n, and such that /(I) < • • • < /(n), f{n -|- 1 ) < • • • < f{n + m). 
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Let US write now T for the tensor algebra over the algebra S of semimartingales, 
equipped with the * product defined in ([71) , so that from now on denotes a 

tensor product of semimartingales in 5®” C T, and J Xi ■ ■ ■ Xn the corresponding 
iterated stochastic integral. We finally obtain the analog for iterated integrals of 
semimartingales of the usual Chen formulas of iterated integrals. Recall that the 
latter hold for either Stratonovich or indefinite Riemann integrals. 


Proposition 4.1. The product of two iterated stochastic integrals is given by: 



where, as above, / runs over all surjections from the set [n + m] to the set [k], for 
max(n, m) < k < m+n, and such that /(I) < • • • < /(n), /(n+1) < • • • < /(n+m). 


For example, the product of two iterated stochastic integrals gives 

J X,X2-J ^3 = J {X,X2X3+X,X3X2+X3X,X2 + {Xi*X3)X2+Xi{X2*X3)). 

The proposition follows from the observation that the inductive rules for the 
quasi-shuffle product in the tensor algebra give the pattern obeyed by products 
of iterated integrals of semimartingales. Namely, setting for X := f A-dB and 
¥■■=/ C-dD, we have 




J {A{ J C_dD))_dB, 

A_dB)C)_dD, 


X»Y 


{AC)-d[B,D]. 


5 . Surjections 

This section presents a concise and mostly self-contained account on the modern 
algebraic theory of surjections, originating independently from F. Hivert’s Ph.D.- 
thesis [TS] and from the work by Chapoton on the permutohedron [3]. We refer to 
the works [151 [24] for more details on the subject. 

Let us write Sj„ ^ for the set of surjective maps from the set [n] n} to 

the set [p], Sj„ p for its linear span, Sj„ for the union of the Sj„ p, p < n, and 

l<p<n 

The linear span of all surjections is denoted 

Sj:=0Sjn,P- 

n,p 

and will be called from now on the set of surjective functions. 

Let us mention for completeness sake that surjective functions in this sense are 
often referred to as word quasisymmetric functions in the literature on algebraic 
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combinatoric. This is because they can be encoded by formal sums of words over 
an ordered alphabet. The set Sj is then written WQSym (this is the notation 
used in the articles we quoted for further details on the underying theory). This 
interpretation, which we will not use in this work, permits to deduce automatically 
certain properties for Sj from general properties of words. See the references [241IT5] 
for more details. 

The vector space Sj is naturally equipped with a Hopf algebra structure, through 
its action on quasi-shuffle algebras |24] . However, we will make use here only of 
the algebra structure. As we just saw, quasi-shuffle algebras are closely related 
to stochastic integration, and this is the reason why Sj will prove to provide an 
appropriate algebraic framework for what is going to be presented in the next 
sections. 

Let us consider a word w over the integers (a sequence of integers) whose set 
of letters is / = {H, •■•,*«} (e.g. w = 35731, / = {1,3,5,7}). Let us write / 
for the unique increasing bijection that sends / = {ii,... ,*„} to { 1 ,... ,n} (e.g. 
/(I) = 1, /(3) = 2, /(5) = 3, /(7) = 4). The packing map pack is the induced map 
on words (for example pack(35731) := /(3)/(5)/(7)/(3)/(l) = 23421). 

Definition 5.1. The product in Sj of f G Sj„ with g G Sj^ ; is defined hy fog := 
h, where h runs over the elements in max(/, k) < i + j < I + k such 

that 

pack(/i(l) • • • h{n)) = /(I) • ■ • fin), pack(/i(n -|- 1) • • • /i(n -I- m)) = 5 ( 1 ) • • ■ g{m). 

This product is associative and unital (the unit can be understood as the unique 
surjection from the emptyset 0 =: [0] to itself), but it is not commutative. 

Proposition 5.2. The linear span of all surjection, Sj, eguipped with the o product, 
is an associative, unital, non-commutative algebra. 

Associativity follows by noticing that the product f o go j of three surjections, 
where j G Sj^ is obtained as the sum of all surjections h with pacfc(/i(l) • ■ • hin)) = 
/(I) • • • fin), packihin+1) ■ • ■ hin+m)) = g(l) • • • gim), packihin+m+V) ■ • ■ /i(n-|- 
m+p)) = j(l) • •-Jb)- 

6. Descents and a Ito-type BCH formula 

It turns out that, similar to the classical theory of iterated integrals, the most 
interesting computations that will take place later in this article on iterated sto¬ 
chastic integrals do not involve the full algebra Sj, but only a small subalgebra, 
known as the descent algebra or algebra of noncommutative symmetric functions 
Sym. For more details the reader is refereed to the standard references [iTl [27] . 

As an algebra, Sym is the free graded associative unital algebra over generators 
In, indexed by non-negative integers. We denote the product * and set ln,m ■= 
In * Im- In general, for a sequence n := ni,... ,nk of integers, we write In '■= 
l„i • • • = 1 = In^. ■ As a vector space, Sym is simply the linear span of the In- 

A surjection / in Sj„ is said to have a descent in position f < n if and only if 
/(*) > /(* + !)• The set of all descents of / is written Desc(/) and 

Desc(/) := {* < n, /(i) > /(I -b 1)}. 

The number dif) that appears in Theorem 1 1.1 1 is the number of descents of /. 
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We also set, for I C [n — 1], 

Desc 7 := 
■= 

and 

Descg/ := 


{/ e Sj„, Desc{/) = /}, 

E /eSj, 

/^Descj 


{/ e Sj„, Desc(/) C /}, 

E /esj. 

/eDescc/ 


When the value of n is obvious from the context, we will abbreviate D" by £>/, 
and similarly for other symbols. Notice that Descc/ = Sjc/Descj, so that, by 
Mobius inversion in the poset of subsets of the set [n — 1], 

(13) Desc/ = y^(—l)hl-NlDesccj. 

JCI 


The subsets Desc/ form a decomposition of Sj„ into a family of disjoint subsets, 
from which it follows that Dj and (by a triangularity argument) D<zi form two 
linearly independent families in Sj„. 


Lemma 6.1. The map i from Sym to Sj defined by: 

^(In) ■ ^C{ni,r!,i+n2,...,r!,iH-|-nfe_i} 

is an injective algebra map from Sym into Sj. 


Proof. Injectivity follows immediately from the linear independency of the Dci- Let 
us show that t is an algebra map. For arbitrary ri = ni,... ,nk, Wi = mi, ... ,mi, 
we have 


'-(In * Im) 


D 


C{ni,ni+n2,...,niH-hrifc.niH-hrafc+mi,....niH-hrifc+miH-hm-i-i} ■ 


On the other hand, L{ln) is, by definition of the o product in Sj, the sum of 

all surjections / G Sj„^+...+„,,such that pack(/(l) ■ ■■ fljii H - h Uk)) 

lies in 


DeSCc{ni,rai+ra 2 ,..,™iH-hrifc-i} 

and pack(/(ni + • • • + + 1) • • • /(ni + • • • + rife + mi + • • • + mi)) lies in 


DeSCc{mi ,mi+m2..,miH-—i}' 

Since there is no constraint on the relative values of /(ni + • • ■ + rife) and /(rii + 
• ■ • + rife + 1), the statement of the lemma follows. g 


The following theorem is the equivalent, in the quasi-shuffle framework, of the 
classical continuous Baker-Campbell-Hausdorff theorem, which computes, among 
others, the logarithm of the solution of a - matrix-valued - linear differential equa¬ 
tion. Theorem 16.21 will appear to play the same role for matrix stochastic linear 
differential equations. It was first stated in [53], but the proof given in that article 
is indirect and relies on structure arguments from the theory of noncommutative 
symmetric functions. Stating those results, which are scattered in the literature 
on algebraic combinatorics, would go beyond the scope of this work. We propose 
therefore a simple and self-contained proof, which is reminiscent of the solution to 
the classical Baker-Campbell-Hausdorff problem stated in m- 
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Let US set I := =• where we write pn for the identity map 

of the set [n] viewed as an element of Sj„. 

Theorem 6.2. We have, in Sj, 


log(-f) = I: i: 


n=l /C[n-1] 


|/| + 1 


C/ 


= E E 


n=l /C[n-1] 


(—1)1^1 fn — 1 

\I\ 


D^. 


Proof. The first part of the statement follows from the computation of the logarithm 
of SyrH) Emd from the previous Lemma. Indeed 

QQ QQ / i\n—1 ^ _ / i\fc—1 

'«e(EL)=E^(i:i.r= E 

n—0 i—1 i—1 

Let us now expand the second term of the identity in the theorem in terms of 
the Df. This yields the coefficient of Df, I = {ii,... ,ik} given by: 

n—k—1 


^ _ "'A" (-ly^ fn-k-l 

^ I TI 4- 1 “ ^ 


Ul + l ^ 7 + fc+l 

ICJGln-l] ' ' j=0 •' 


n — k — 1 


n — k — 1 

j 




= (- 1 )'^ E 

1=0 { 


from which the theorem follows. 


In view of Proposition 14. II Theorem 1 1.1 1 also follows. 

The elements E/c[n-i] ' (”s|^)” ^ Df are analogs in Sj of the celebrated 

Solomon idempotents, see [23] for further details. 

7 . NONCOMMUTATIVE stochastic CALCULUS 

In the previous sections of the article, we investigated quasi-shuffie-type proper¬ 
ties of iterated integrals of semimartingales. We also showed the strong relationship 
between them and properties of surjections. Now Proposition 14.II is restated and 
generalized: 

Proposition 7.1. The product of k iterated stochastic integrals of semimartingales 
is given by: 

. 

/eDescc{„i,...,„i + ...+„j^_j} '' 
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More generally, for fi S ,... 



/i( J XI---XI' 

)■ --- ■fk{J x^- 

\rk 

-^nk 

0 

0 

II 

fk){ J 

--x^- 


In the last formula, from which the first follows, it is implicitly assumed that the 
action of Sj^. on f Xi ■ ■ ■ Xk is extended linearly to the linear span of Sj^,, that is, 
the action of a linear combination of surjections fi is the linear combination of the 
actions of the fi. 

We let the reader check that the k = 2 case of the last formula follows from 
the definition of the f(^J Xl ■ ■ ■ X^'^ and from Proposition 14.11 The general case 
follows by induction. 

For example, the expansion of the triple product (/ X) ■ (f Y 1 Y 2 ) ■ (f Z 1 Z 2 ) 
includes terms such as: 


J XYi{Y 2 *Zi)Z 2 , j {X i.Yi){Y 2 i^ Zi)Z 2 and J Zi(X * W * ^ 2 )^ 2 - 


The purpose of the present section is to extend this picture to the operator 
setting, that is, to iterated stochastic integrals of, say, n x n square matrices M = 
M{X^’^)i<ij<n whose entries are semimartingales. Note that we write 

the indices of the entries as exponents for notational convenience in forthcoming 
computations. The set of such matrices is denoted 971. 

For matrices Mi ,..., Mk £ 971, we set: 

/ Ml ■ • • Mfe := ( V [ M*’*®M*®’*" • ■ • 

Zi 

and for / e Sj^, 

fH Mi---Mk) := ( J 


Last, if Xifi is a linear combination of surjections in Sj;,, then we set: 

j Ml---Mk) :=^A,/,( J Ml---Mk). 


The product rule of Proposition 17.11 applies (entry-wise), and we obtain, for 
Ml ,..., Mk ,Ni,..., Nra S 971: 

f f N,..N„ 


= ( E , E 

/eDesCc{fc} + m 


f{M, 




)*fc + l 5^fe4-2 


= f{J Mi---MkNi---NAj. 

/GDesccifc) 


m 
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For example, let fc = 2,m = 1, and consider 2x2 matrices, this gives for the 
first entry of the product: 

( J Ml M2 • J ^ ^ 


L 1<2 


+ * N('^) + (M^* TVf . 

For higher products we obtain similarly: 

Proposition 7.2. For Ml ,..., M^^, Mf,..., M^^ G dJl, we have: 

( J Ml---Ml^) .(y M^.-.M„\) 


E 


/( I Mi.-.Mi^..-M^.-.M„\) 


/eDesCc{ni,...,„i + ... + nj^_j} 

and more generally, for fi G Sj„ fk G Sj„ we obtain: 


/i(y Mi-.-Mi^) •... •/fe(y m^...m„\) 

= (/iO---o/fe)( /M1..-M1^..-M^--.M„\). 


The proposition follows from the linear case (Proposition 17.ip by expanding 
entry-wise the products of matrices. 

We are now in the position to calculate the logarithm of the ltd-Taylor series. 

Theorem 7.3. For an arbitrary matrix M G 9Jt and X = I (formal) 

solution of the stochastic differential equation dX = X-dM, Xq := Id, we have: 


\og{x) = Y, E 

n=l /C[n-1] 


(- 1 ) 1^1 


n — 1 

|/| 




M”. 


This formula may provide the basis for interesting numerical properties, 
instance, truncating the expansion of log(X) at order fc, that is, looking at 
k / 1 \in / i\ -1 


For 


E E 

n=l /C[n-1] 


(- 1 ) 1^1 


n — 1 

1^1 


D) 


M" 


and applying the exponential map, can be expected (in view of similar phenomena in 
the deterministic case) to provide a better approximation to X, than the truncation 
of the original expansion f 

The first few terms of the expansion of log(X) read: 

log(X) = Jmf (^^12) - ^(21) + (11))^ J m2 

+ (i(123) - i((213) + (312) + (112) + (212)) 

-i((132) + (231) + (122) + (121)) 

^(321) (211) -h (111) (221))) J + - . 
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We remind the reader that we represent a surjection / G Sj^. by the sequence of its 
values (/(l)---/(fc)). 

For example, for 2x2 matrices, the (l,2)-entry of the term (11) f M'^ is given 
by E^<2f the one of (231) f reads E^,J<2 f the 

one of'(221)/M3 reads E*.j<2 / 
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